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Abstract 

We introduce and investigate the wellposedness of two models describing the self- 
propelled motion of a “small bio-mimetic swimmer” in the 2-D and 3-D incompressible 
fluids modeled by the Navier-Stokes equations. It is assumed that the swimmer’s body 
consists of finitely many subsequently connected parts, identified with the fluid they oc¬ 
cupy, linked by the rotational and elastic forces. The swimmer employs the change of its 
shape, inflicted by respective explicit internal forces, as the means for self-propulsion in a 
surrounding medium. Similar models were previously investigated in [15]-[19] where the 
fluid was modeled by the liner nonstationary Stokes equations. Such models are of interest 
in biological and engineering applications dealing with the study and design of propulsion 
systems in fluids and air. 
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1 Introduction 

The swimming phenomenon has been the subject of interest for many researchers in various areas 
of natural sciences for a long time, aimed primarily at understanding biomechanics of swimming 
locomotion of biological organisms, see Gray [T^(1932), Gray and Hancock [TS] (1951), Taylor 
pH] (1951), [35] (1952), Wu [33] (1971), Lighthill [33] (1975), and others. This research resulted 
in the derivation of a number of mathematical models for swimming motion in the (whole) or 
M^-spaces with the swimmer to be used as the reference frame, see, e.g., Childress [5] (1981) and 
the references therein. In particular, based on the size of Reynolds number, it was suggested (for 
the purpose of simplification) to divide swimming models into three groups: microswimmers 
(such as flagella, spermatozoa, etc.) with “insignificant” inertia; “regular” swimmers (fish, 
dolphins, humans, etc.), whose motion takes into account both viscosity of fluid and inertia; 
and Euler’s swimmers, in which case viscosity is to be "neglected”. 

‘This work was supported in part by the INdAM national group GNAMPA and the INdAM-CNRS GDRE 
CONEDP, and by the Grant 317297 from Simons Foundation. 

^The work by this author was supported in part by NSF Grant DMS-1007981. 
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It also appears that the following two, in fact, mutually excluding approaches were distinguished 
to model the swimming phenomenon (see, e.g., Childress, [S] (1981)). One, which we can call the 
shape-transformation approach” , exploits the idea that the swimmer’s shape transformations 
during the actual swimming process can be viewed as a set-valued map in time (see the seminal 
paper by Shapere and Wilczeck [5^ (1989)). The respective models describe the swimmer’s 
position in a fluid via the aforementioned maps, see, e.g., |5] (1981), [53] (2008), [7] (2011) 
and the references therein. Typically, such models consider these maps as a priori prescribed, 
in which case the question whether the respective maps are admissible, i.e., compatible with 
the principle of self-propulsion of swimming locomotion or not, remains unanswered. In other 
words, one cannot guarantee that the model at hand describes the respective motion as a self- 
propulsive, i.e., swimming process. To ensure the positive answer to this question one needs to 
be able to answer the question whether the a priori prescribed body changes of swimmer’s shape 
can indeed be a result of actions of its internal forces under unknown in advance interaction 
with the surrounding medium. 

The other modeling approach (we will call it the “swimmer’s internal forces approach” or SIF 
approach) assumes that the available internal swimmer’s forces are explicitly described in the 
model equations and, therefore, they determine the resulting swimming motion. In particular, 
these forces will define the respective swimmer’s shape transformations in time as a result of 
an unknown-in-advance interaction of swimmer’s body with the surrounding medium under the 
action of the aforementioned forces. For this approach, we refer to Peskin [23] (1975), Fauci 
and Peskin |5] (1988), Fauci jH] (1993), Peskin and McQueen [21| (1994), Tytell, Fauci et al [JT] 
(2010), Khapalov [T7| and the references therein. 

The original idea of Peskin’s approach is to view a “narrow” swimmer as an immaterial “im¬ 
mersed boundary”. Within this approach the swimming motion is defined at each moment 
of time by the explicit swimmer’s internal forces. Following the ideas of Peskin’s approach, 
Khapalov introduced the immersed body SIF modeling approach in which the bodies of “small” 
flexible swimmers are assumed to be identified with the fluid within their shapes, see m-m 
(2005-2014). Indeed, in the framework of Peskin’s method the swimmer is modeled as an im¬ 
material curve, identified with the fluid, further discretized for computational purposes on some 
grid as a collection of finitely many “cells”, which in turn can be viewed as an immerse body, 
see Figures 1 and 2. The idea here is to try, making use of mathematical simplifications of 
such approach, to focus on the issue of macro dynamics of a swimmer. The simplifications 
(they seem to us to be legitimate within the framework of our interest) include the reduction 
of the number of model equations and avoiding the analysis of micro level interaction between 
a “solid” swimmer’s body surface and fluid. It should be noted along these lines that in typical 
swimming models dealing with “solid” swimmers, the latter are modeled as “traveling holes” in 
system’s space domain, that is, the aforementioned “micro level” surface interaction is not in 
the picture as well. 

In the above-cited works by Khapalov m-m, the immerse body SIF approach was applied to 
the nonstationary Stokes equations in 2-D and 3-D dimensions with the goal to investigate the 
well-posedness of respective models and their controllability properties. In this paper our goal 
is to extend these results with respect to well-posedness to the case of Navier-Stokes equations 
in both the 2-D and 3-D incompressible fluids. To our knowledge, there were no previous 
publications investigating this issue within the SIF approach. 

Related references on well-posedness of swimming models. To our knowledge, in the 
context of PDF approach to swimming modeling, the classical mathematical issues of well- 
posedness were addressed for the first time by Galdi [10] (1999), [H] (2002) for a model of 
swimming micromotions in (with the swimmer as the reference frame). In |25j (2008) San 
Martin et al discussed the well-posedness of a 2-D swimming model within the framework of 
the shape transformation approach for the fluid governed by the 2-D Navier-Stokes equations. 
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Figure 1: 2-D swimmer consisting of 4 rectangles. 

Swimming models in the framework of ODE’s. A number of attempts were made to 
introduce various reduction techniques to convert swimming model equations into systems of 
ode’s (e.g., by making use of empiric observations and experimental data, etc.), see, e.g., 
Becker et al [3] (2003); Kanso et al. [Hj (2005); Alouges et al. [T] (2008), Dal Maso et al. [7] 
(2011) and the references therein. 

The paper is organized as follows. In Section 2 we state our main results. In Section 3 we 
discuss in detail the modeling approach of this paper to swimming locomotion. In Section 4 we 
prove our main results after stating several auxiliary lemmas, proven further in Section 5. In 
Sections 6 and 7 (Appendices A and B) we remind the reader some classical results used in our 
proofs. 


2 Main results 

To formulate these results, we will need the following function spaces. 

Let D C be a bounded domain with locally Lipschitz boundary dfl. Below, we use the 
following classical notations: 
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Figure 2: 3-D swimmer consisting of 4 parallelepipeds. 


• d denotes the dimension of the space domain, equal either to 2 or to 3; 

• C“(D) denotes the space of infinitely many times differentiable functions with compact 
support in fl; 

• V'{^1) denotes space of distributions in il, i.e., the dual space of C“(f2); 

• VF^’P(r2), 1 < p < oo denote the Sobolev spaces over D, i.e., the Banach spaces of functions 
in LP(D) whose first (generalized) derivatives belong to LP(D); 

• H^{n) = and H^{n) = 

• Ffo(D) denotes the subspace of consisting of functions vanishing on dfi. 

denotes the dual space of 

Following m, page 5, we also introduce the following d-dimensional vector function spaces: 

V:={pG ; divp = 0}, 

H := c1l2(V) , V := cl^i(V) = {p G [ido(D)]‘' ; divp = 0} , 
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where the symbol c1l 2 stands for the closure with respect to the [L^(n)]‘^-norm, and cl/^i - with 
respect to the [i?Q (n)]'^-norm. The latter is induced by the scalar product 

d dp 

<(/?,'(/'> := ^ hpixjtpix^dx, (p= ')p = 

i=i i j=i o 

where Dj is the differentiation operator with respect to Xj. To simplify notations, below we will 
use the notation ||<p||l 2 (resp. ||<p||//i) both for functions (p S L^(n) (resp. p S Hq{V,)) and for 
functions p e [L^(n)]'^ (resp. p e [iJQ(n)]‘^). 

Let V and H' stand for the dual spaces respectively of V and iJ. Then, identifying JI with 
H' , we have 

V C H = H' cV . 


Our main results deal with the well-posedness of 2-D and 3-D swimming models (§, @- 

(14), described in detail in the next section and visualized by Figures 1 and 2. 


Theorem 1 (Well-posedness of the 2-D swimming model) Letd = 2 and for some T > 0 
let uq G H, • I znp G Ki,..., kn-Ii Vi, , vn -2 G L^(0, T). Assume that the assump¬ 

tions (H1)-(H2) (given in Section 3) hold, and that 

\zip - Zjp\ > 2r , i,j = i^j, (1) 


where r > 0 is the constant in (HI). (Assumption (1) ensures that no parts of swimmer’s 
body overlap with each other and all lie within VL.) Then, there exists T* G (0,T] such that sys¬ 
tem (PI), (p]), (T^,(ll )/{12),(13) admits a unique solution (u, z) in C(0, T* ; iL)nL^(0, T* ;V)x 
[C([0,r*];M‘^)P, and 


S(zi{t))cn, \zi{t) - Zj{t)\ > 2r, i,j = i^j VtG[0,T*]. (2) 

The formula for Vp, complementing the given pair {u, z) in models is given 

in Proposition\6.1\ below. 


Theorem 2 (Well-posedness of the 3-D swimming model) Let d = 3, dLl be of class 
and Uq gV . Then the result stated in Theorem 1 holds for model 0), pp, @ 

for a unique triplet (u,p,z) such that ut,Au,'\/p G [L^(Qt*)]^ and u G (^([O, T*]; F), where 
Qt. = (0,T*) X O. 


Theorem 3 (Additonal regularity) Let dLl be of class . If uq G [iL^(fl)]^ f) F, then in 
Theorem^solution u lies in [iL^’^(Qr.)]^ P|C'([0,T*]; F). In turn, if uq G [iL^(O)]^ P| F, then 
uG [H‘^’^{QT>)f C]C{[Q,T*]-V) in Theorem\^ 

Here, = {p \ p G L^(0, T*; iJ^(O)), pt G L^(Qt*)}. This result is an immediate 

consequence of Theorem in Section 6 (Appendix A). 

Remark 2.1 The duration of time T* (i.e., of existence of solutions) in Theorems^Tj^ depends 
on the parameters of respective model uq, Zip’s, ki’s, Vi’s and the initial shape and position of 
the swimmer in Ll. One can view T* as a new initial instant of time to apply these theorems 
again to further extend the interval of existence from T* to some T** > T*, provided that 
the assumptions of the corresponding theorem hold at t = T*, and so on. Our proofs below 
indicate that in the 2-D case one can extend this time up to the collision of swimmer with the 
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boundary of Q or up to the moment when swimmer’s body parts will collide with each other, 
see conditions In turn, these circumstances depend on or can be regulated by a suitable 

choice of swimmer’s internal forces, i.e., functions Vi’s (and ki’s when they can vary, see the 
next section). To the contrary, in the 3-D case, such interval of existence of solution will also 
depend on uq, Zi^ ’s, ki’s, and Vi’s via an additional condition Theorem^ (compare it to 

Theorem^. 


3 Swimming model 

Following HZI, ng-ni], we describe the locomotion of a swimmer in a fluid by a hybrid nonlinear 
system of two sets of pde/ode equations: 

{ Ut — i^Au + {u ■ V)u + Vp = / 
div u = 0 
u = 0 

u(0,-) = Uq 


in (0,r) X O, 
in (0,T) X n, 

in (o,r) X an, 

in fl, 


dzi 

dt 


1 

meas(S'(0)) 


/S(z.(t)) 


u{t, x) dx, 


Zi{t)) = Zifi, 


i = te(0,T). (4) 


System ^ describes the evolution of incompressible fluid due to the Navier-Stokes equations 
under the influence of the forcing term f{t, x) representing the actions of swimmer. Here, 
(u{t,x),p(t,x)) are respectively the velocity of the fluid and its pressure at point x at time 
t, and v is the kinematic viscosity constant. In turn, system describes the motion of the 
swimmer in H, whose flexible body consists of N sets S{zi{t)) within H. These sets are identified 
with the fluid within the space they occupy at time t and are linked between themselves by the 
rotational and elastic forces as illustrated on Figures 1 and 2. The points Zi{t)’s represent the 
centers of mass of the respective parts of swimmer’s body. The instantaneous velocity of each 
part is calculated as the average fluid velocity within it at time t. 


We will now describe the assumptions on the parameters of model @-0 in detail. 


3.1 Swimmer’s body 

Below, for simplicity of notations, we will denote the sets S{zi{f)),i = 1,... ,N also as S{zi) or 
Si{zi). Throughout the paper we assume the following two main assumptions: 

(HI) All sets S{zi),i = 1,... ,7V are obtained by shifting the same set S'(O) C i.e., 

S{zi) = ZiS{0), i = l,...,N, (5) 

where S{0) is open and lies in a ball Br{0) of radius r > 0, and its center of mass is the 
origin. 

Remark 3.1 The results of this paper will hold at no extra cost if we assume that the swimmer 
at hand consists of different body parts Si{0) C Br{0), in which case one will need to replace 
0 with S{zi) = Zi + 5'i(0), i = 1,...,N and add respective normalizing coefficients in the 
expressions for swimmer’s internal forces to ensure that they satisfy the 3rd Newton’s 

Law. In particular, the swimmers on Figures 1 and 2 consist of identical sets each of which has 
its own orientation in space. 
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We will need the following concepts from Section 3.11] to formulate our second assumption 
on the geometry of swimmer in this paper, which will also be used in the respective proofs 
below. 

Given a bounded set fl G and a unit vector v G M", denote by 

TTy G K" ; ^ = 0} , flj, := {^ G TTy ; 31 G K s.t. ^ +tv € ft} 

respectively the orthogonal hyperplane to v and the orthogonal projection of fl on tt^. Then, 
for every y G we will call the set 

:= {t G M ; y + tv G Q, y G fli,} ( 6 ) 

the section of fl corresponding to y. Accordingly, given a function (/i: 17 —> K and any y G 17^, 
we define the function (^^^{t), called section of (p corresponding to y, as 

4>l-.n^9y^t^4>i{t)&R, <t>y{t)-.= ct>{y + tv) . (7) 

(H2) There exist positive constants Hq and /Cg such that for any vector h G \ {0} we 

can find a vector rj = r]{hPj, | 77 |= 1 which satisfies 

meas(S'A)H = [ dt < ]Cs\h\ V j/ G 17^ , ( 8 ) 

7(Sa)? 

where Sa ■= {h + S(0)) A S(0) is the symmetric difference between S'(O) and h + S{0), i.e. 
Sa = {{h + SiO)) U 3(0)) \ {{h + 5(0)) n 5(0)). 

Assumption (H2) means that size of the projection of 5 a on the hyperplane in R'^ perpendicular 
to vector 77 changes uniformly Lipschitz continuously relative to the magnitude of the shift h 
of the set 5(0) in the direction of h. This assumption is principally weaker than the respective 
assumption on the regularity of the shifts of 5(0) in [T7], [IS]-[Tn], where 77 was always selected 
to be h (we will illustrate it in Example 3.1 below). In the case when 77 = h, it is easy to verify 
that (H2) is satisfied for discs and rectangles in 2-D and for balls and parallelepipeds in 3-17. 


Remark 3.2 In this paper we assume that all swimmer’s body parts are identical sets. One 
can choose these sets to be of distinct shapes and sizes, in which case, however, the respective 
normalizing coefficients should be added to the forcing terms to ensure that all swimmer’s forces 
are to be its internal forces. 


We conclude this subsection with an example showing that there exist particular shapes of the 
set 5(0) which require the presence of 77 = 77 (h) in (H2) instead of the straightforward choice 



Example 3.1. Fix a constant k > 0. We claim that there exist sets 5(0) which satisfies (H2) 
for some 77 , but 77 cannot be selected to be co-linear with h. 

Let us consider the set 5 shown on Figure 3, 


S ■■= {(a;, 77 ) G I X G (0,1), 7 / G (a(x),/3(x)), } 
where the functions a, (3: [ 0 , 1 ] —>■ K are defined as follows: 

f —X, if 0 < X < 
a(x) := < 


1 

2 ’ 


/3(x) := < 


2k X — 


2"+i 


if X G 


X — 1, 


if 1 > X > 


1 

2 ’ 



if X G 


1 3 

2n+l ’ 2'^+2 ’ 

3 1 ■ 

2Tt + 2 ’ 2’T’ ’ 


7 













2" + l 



Figure 3: Set S in Example 3.1. 


n = 0,1,.... Observe that the part of the boundary dS corresponding to y = /3{x) is given 
by lines with slopes either 2k or —2k. If we introduce the notation Tm {m > 1) for the 
isosceles triangle of base 2“"* in [2“"*, 2“™+^] and height k/2 ’”, and To for the isosceles triangle 
{(a;, y) G ; a: G [0,1], y G [Q!(a:), 0]}, then we have 

5 = \jTm- 

771>0 

Denote by b{S) the center of mass of S and set S{0) = S — b(S'). We claim that, no matter 
what ho and JCs we choose, if we shift 5'(0) hy h = £(—1,0), for a suitable £ G (0, ho), and if we 
use ri{h) = (—1,0), then we can always find y G D,, such that 

meas [S{0) A {h + 5'(0))]^ > ICs e, 

and, thus, ® does not hold. On the other hand, by setting r][h) = (0,1) for all h G B{0, 1), we 
can prove that 

meas [S'(O) A {h + S'(O))]^ < 4k£ V y G D,, , (9) 

so that (|^ is satisfied with ho = I and JCs = 4k. 

Since the Lebesgue measure is invariant with respect to translations, in the computations below 
we will always use S and Sa = S A (h + S') in place of S(0) and S(0) A {h + S(0)). 

We start with the negative result. For every fixed ho > 0 and JCs > 0, let us consider m large 
enough to have 2“”^ < ho and 2(m — 1) > JCs- By choosing h = 2“’”(—1,0) G Bho{0), we 
claim we can find y G (Sa)^ C D/j such that meas(SA)^ > JCs \h\ = Ksl2J^, so that the choice 
y(h) = is not suited for this set S. 

The projection (Sa)/! on the y-axis (which is the orthogonal line to h) is the segment (—1/2, k/ 2). 
Thus, the value y = 2 “™k belongs to (Sa)/i and (see Figures 3- 4) 

{.Sa)1^ U 

l</c<m 
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i.e., (Sa)^ consists of the union of the sections of the symmetric differences of the triangles in 
S whose height exceeds 2~'^k. Hence, meas (S'a)^} >= 2(m — l)|/i| > K-s \h\ from which our 
claim follows. 



Figure 4: The choice of section line in Example 3.1. 


We now show that for e < 1, /i = £(—1,0) and 77 = ? 7 ((—1,0)) = (0,1), the sections of the 
symmetric difference Sa for all y G satisfy ([^. 

Let h = £(— 1,0). Consider, as a preliminary step, the symmetric difference 7a = %n A (h+Tm) 
for any fixed m S N. It is not difficult to verify that for y G ftrj 

meas(7A)^ < 2 «;£ 

as illustrated by Figure 5. To pass to the sections of Sa, it is sufficient to observe that (<5'A)r; 
is the interval (—£, 1) on the x-axis and that or any y G (Sa)?] C fl?/, there holds: 

(5A)^C(roA(h + ro))^u(r^A(h + r^))^, when ye ( 0 , 1 ), 

('S'a)^ A {(a, &) e ; a e [—£, 0 ] , 6 e [—2tt;a, 2 ko]}^ , when y e (—£, 0 ]. 


3.2 Swimmer’s internal forces 

In this subsection we give the precise description of the force term / in ([^. Due to the nature of 
the swimming motion as self-propulsion, all the forces in model are internal relative to the 

swimmer, i.e., their sum is equal to zero and their torque is constant. In turn, these forces, acting 
between swimmer’s body parts, will create a pressure upon the surrounding fluid and, thus, will 
act as external forces upon it. We assume that all forces act through the immaterial links 
attached to the centers of mass Zi{tys of the sets S(zi{t)), and then are uniformly transmitted 
to all points in their respective supports. 

Similar to m, in this paper we consider two types of forces forming / in ([^: rotational forces 
and elastic forces, which we represent as 

f{t,x) := frot{t,x) + fel{t,x) . (10) 
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e 



Figure 5: The measure of the sections (7 a)^ can be estimated by the height of the right triangle 
with base £ and hypotenuse parallel to the side of Tm- Left: e > Right: e < 


The 2nd term in (10) describes the forces acting as elastic links between any two subsequent 
sets S'(zi)’s to preserve the integrity of swimmer’s structure. They act according to the 3rd 
Newton’s law and Hooke’s law with variable (positive) rigidity coefficients Ki(t),..., KAr_i(t) 
when the distances between any two adjacent points Zi{t) and i = 2,...,N, deviate 

from the respective given values it > 0,i = 1,..., N — l(see [TTjl: 


N r 


fel{t,x) := 


i- U \ tj.\ •2i-l(t)| (-i-l ( f ^ f 

-- (zRf) - Zi-i[t)) 


■ ^iit,x) Ki_i{t) 


\Zi{t) - Zi-i{t)\ 

\zi{t) -Zi_i(t)| -£*_1 


\Zi{t) - Zi_l(t)| 

Here ^i{t,x) denotes the characteristic function of S{zi{t)), i.e., 


[zi-iit) - z^{t)) 


( 11 ) 


^tit,x) := 


1 if a; S S{zi{t)) 
0 otherwise . 


Remark 3.3 In the above structure we can also assume that Ki(t),... ,KN-i{t) can be of any 
sign and replace (11) with more general (and simpler) description of pairs of co-linear forces 
between Zi ’s; 


AT r 


fei{t,x) := Y 


fi-l{t,x) K^-l{t) (zi(t) - Zj_i(t)) 
fiit,x) {zi-i{t) - Z^(t)) 


( 12 ) 


All the proofs in this paper are given for technically more complex case of elastic forces in (11). 


The 1st term in (10) describes the forces that allow each point Zi{t),i = 2,..., N— 1 to rotate any 
pair of the adjacent points Zi-i(t} and Zi-^.l(t) about it in either folding or unfolding fashion. In 
turn, by the 3rd Newton’s law, these points will act back on Zi(t) with the respective countering 
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force. The description of rotational forces requires principally different approaches for the 2-D 
and 3-D cases. 

In the 2-D case all the forces lie in the same plane at all times, and we can describe them 
by making use of the matrix 


as follows (see my- 



N-l 


frot,2di*,x) ■■= ^ 




'Zi-i{t) - Zi{t)\'^ 


\Zi+i(t) - 


A (zi+i(t) - Zi{t)) 


N-l 


+ X! Vi-lit) 


A [zi{t) - Zi-iit)) 


\Zi-i{t) - Zi{t)\^ 

\z^+lit) - 


A (zj(f) 



(13) 


where functions Ui(t),... ,ujv-2(0 characterize the strength and orientation (folding or unfold¬ 
ing) of respective pairs of rotational forces at time t. 

In turn, in the 3-D case , to satisfy the 3rd Newton’s law, we need to make sure that the 
respective rotational forces acting on Zi-iit) and Zi+iit) lie in the same plane spanned by 
the vectors Zi_i(t) — Zi{t) and Zi+i(t) — Ziit). In order to achieve the continuity of these 
forces in time, in this paper we choose to reduce their magnitudes to zero, when the triplet 
{zi-i{t), Zi{t), Zi+i{t)} approaches the aligned configuration (for other options see [IH])- Indeed, 
such configuration admits infinitely many planes containing this triplet, which makes it an 
intrinsic point of discontinuity for the procedure of the choice of the rotational plane by means 
of the rotational forces whose magnitudes are strictly separated from zero. Respectively, we 
define the 3d rotational forces as follows: 


N-l 


/rot,3d(t,a;) := X! D-i{t) 




- it+iit,x) 


^i-lit,x) Pi[t]{z^-i{t) - Ziit)) 
'z^-iit) - Zi(f)P 


\Zi+iit) - Zi(t)|" 


Qi[t\{zi+iit) - Ziit)) 


N-l 


+ X! Vi-lit) 


Pi[t]{ziit) - Zi-lit)) 


\zi-iit) - Ziit)y 


Q^[t]{ziit) - Zi+lit)) 


(14) 


\Zi+iit) - 2:j(f)|2 

where the scalar functions Ui(t),..., ujv- 2(0 control the magnitudes of the rotational forces and 
determine whether they act in folding or unfolding fashion, and 


X Pi[t]x := [iz,-iit) - Ziit)) X izi+iit) - Ziit))] X X, 

X 1 -^ Qi[t]yi := X X [izi-iit) - Ziit)) x (zi+i(t) - Zi(t))] . 

Note that Pi[t]x = —Qi[t]x and |Ri[^]^l = ^ for any x when points Zi_i(t), Ziit), z i-\-iit) 

converge to the aligned configuration. 
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Remark 3.4 The forcing term / in §». ( |l0| ) can (more precisely) he denoted as f{t, x;z,k,v). 
However, we will use a shorter notation f{t, x) as, typically, there is no ambiguity about the 
choices of z,k and v. 

We can prove the following result. 


Lemma 3.1 Let z = (zi, ..., zn) G [C([0, T] ; C [C([0,T] ; k = (ki, ..., kn-i) G 

[L^(0,T)]^“^ and v = {vi,... ,vn- 2 ) G [L^(0,T)]^“^ be fixed. Assume that for all t G [0,T] 
there holds 


\zi{t) - Zj{t)\ > 2r, 


i,j = i^j, 


with r > 0 as in (HI). Then, the forcing term f{t,x) defined in (|l0| belongs to L^(0, T ; [L^(f2)]‘^) 
and there hold the following estimates 


N 


ll/ei|lL2(o,T;[L2(n)]ti) < 2-^meas(n) || |«i-i|||L^(o,T) 


i=2 


X {ll^i ~ ^i-l|lc([0,T] + ^i-l} 


N-1 


||/roi,2d||L2(0,r;[L2(n)]2) < 2 ■\/meas(ft) || ki-l|HL^(0,T) 


i=2 


X max \\\Zi- 2:^-1 II C([ 0 ,T];R 2 ) + 

2 — Z,..., i V i I 


N-1 


ll/rot, 3d||L2(o,r;[L2(n)]3) < 4 ■\/meas(n) II |^i-i|||l3(o,t) 


X “ 2i-lllc([0.T] ;R3) • 


The above estimates yield that 


l/llL2(0,T;[L2(n)]‘i) < C: 


\Zi - Zi-i\ 


C([0,T] ;R2) 


2 r 


(15) 


where positive constant C > 0 depends on T, 12 and L^(0, T)-norms of parameters Vi’s and Ki’s. 


Remark 3.5 It is not difficult to see that the above-defined forces in (10)-(14) are internal 
relative to the swimmer at hand, i.e., their sum is zero and their angular momentums are 
constant (see, cf. J27[ Chapter 12]). 


4 Proofs of the main results 

The proof of Theorem will be based on Schauder’s fixed point theorem, and will follow from 
a series of lemmas. The general scheme of our proofs is similar to that introduced in m, m, 
[It] (Chapters 11-12), [19] for the case of the fluid governed by the linear nonstationary Stokes 
equations. In this paper we consider the full nonlinear Navier-Stokes equations, which will 
require some principal technical modihcations of this scheme. Our arguments below are nearly 
the same for both 2-D and 3-D cases. However, in the part dealing with the Navier Stokes 
equations, the latter is traditionally more challenging. 
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Remark 4.1 In our our proofs given below we assume that all the parameters in Theorems^ 
and[^ are fixed, namely: T > 0, the initial datum uq either in H (if d = 2) or in V (if d = 3^, 
a vector Zq = (zi.o, • ■ •, ZAr,o) S as in ^ and functions ki, ..., k^-i, vi,...,vn -2 & 

L^(0,r). We may also omit the explicit mentioning of dependence of some of the constants below 
from these parameters. 


We will use the symbol Ti to denote the length of the time-interval of existence of solutions to 
Navier Stokes equations in §■ For d = 2 , due to Theoremcited in Section 6, Ti = -boo. 


For d = 3, the size of Ti > 0 depends on (/, Mq) as cited in Theorem]^ in Section 6 (Appendix 
A). 


More precisely, in this case we will need to identify the value of Ti which, given uq, will work 
uniformly for any selection of / used in the proofs of Theorems I][2 below. Namely, we set: 


Ti G (0, min 


T,C° \\uo\\l^ 



) id = 3), 


(16) 


where ( is from (15). 


In subsection |4.1| we introduce necessary notations and will state the aforementioned lemmas, 
whose proofs will be given in section]^ In subsection 4.2 we will prove Theorems l][2 


4.1 Formulation of auxiliary lemmas 

Set Bq. For any T* G (0,Ti) and q > 0, set 


Bq := 



L^{0,T*-V) 



We will show below that, for suitable choice of pair {T*,q), the following three nonlinear oper¬ 
ators will be well defined. 


Operator T. Set 

T: L2(0,r* ;F) D e, ^ [C([0, T*]; C [C([0, T*]; 
Tu:=z= {zi,...,zn), 


N 


(17) 


where zi, ..., zn are trajectories of system (|^ with the aforementioned fixed initial state zq- 
We will show that for suitable pair (T*, q) they will satisfy ([^ . 


Operator F. For any z G T{Bq), set: 

F: [C([0, T *]; D [C([0, T *]; fl)]^ D T{Bq) - 

Fz = F(zi, ...,zn)-= f, 


L2(0,r*;[L2(fI)f), 


( 18 ) 


where / is the forcing term defined in (10), (13 )-(|14[). 
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Operator S. For any / G FT{Bq), we set respectively: 


S: L2(0,T* ;[L2(^])]‘^) D FT(eg) —^ L2(0,T*;F), 

(19) 

Sf:=u, 

where u is the solution to uncoupled Navier-Stokes equations ([^. 

To prove Theorem 1, we intend to show that the mapping 

SFT:Bg —^ L^{0,T*;V) 

has a unique hxed point in Bq for some pair (T* ,q). To this end, we will show that all the 
above operators are continuous, while T is also compact. Then the desirable result will follow 
by Schauder’s fixed point theorem. 

Lemma 4.1 Fix T > 0, u G L‘^{0,T ;[Hq{Q)]'^), d = 2,3, and assume |^. Then there 
exists T* G (0,T] such that system 0 admits a unique solution z in [C([0, T*]; 11)]^ C 
[C([0,T*];K‘^)]^ and § holds. Finally, given a constant k > 0, there exists G (0, T] such 
that 0) holds on [0,T^] (in place of[Q,T*\)for any choice ofu such that ||u||l 2 (q r.[^i(f 2 )]d) < k. 


Lemma 4.2 Under the assumptions of Lemma \4.1[ for any fixed q> 0, let T* be such that 

21 


0 < T* < min<^ T* 


I meas 


(g(o)) y 


’ V cqCqICs q J 


( 20 ) 


where T* is defined in Lemma 4-1 (see (33)), cq is the constant in ( [51] ), Cq is the constant 
in (52) and ICs is the constant in (H2). Then, operator T defined in 0 is well defined, 
bounded, continuous and compact. In particular, there exist constants L'r,L'j. > 0 such that for 


all u,v G Bq: 


and 


|Tu|| 


[C([0,T*] ;I 


< 


Lt ^I^oI 


iTu- Tul 


Vt*\ 


y/T* 


w|1l2(o,t* -v) I ! 


[C([0.T-] KsqVT^ “ ''^llL2(0.r* -V) ■ 


( 21 ) 

( 22 ) 


1 - 


meas(S(0)) 


Lemma 4.3 For any fixed q > 0 and T* as in Lemma 4-2, operator F m (18) is well defined, 
bounded and continuous. In particular, there exist positive constants Ly,L'^ such that for z,w G 
T{Bq) we have: 

I|F2;||l2(o.t* ;[L2(n)]'i) < Lf r„_„(T*) < (23) 


where the last inequality is due to (15), and 

IjFz — Fx(;||l2(o,t* ;[L2(n)]rf) < LpT,^^y{T*) — wH [c([o,t*] > (24) 

rK;,„(T*) := max |||ft:i||L 2 (o,T*)) • ■ • ; |I^Af-l|lL 2 ( 0 .T*)) lkl||L 2 ( 0 ,T*)) ■ • • j |kAf- 2 ||L 2 ( 0 .T*)} • 


Remark 4.2 The above estimates (23)-(24) can be refined to include explicitly the [C([0,T*]; 
norms of trajectories z,w G T{Bq). However, this is not principal, because T{Bq ) C [C([0,T*] ;n)]^. 
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Lemma 4.4 For any fixed T* € (0,Ti), operator S defined in (191 is well defined, bounded and 
continuous. In particular, there exist constants Ls,Lg and Cs such that (see Remark\4.1]j 


l|S/llL2(o,r* ;y) < Ls (||mo| 


L 2 


IL2(0,T* ;[L2(n)]<i 


)) 


(25) 


and for d = 2: 

IIS/ — Sg||L 2 (Q_-r» .y) 

< L's e^^ .11911=2,0,^. iij _ , (26) 

while for d = 3; 

l|S/ - S5 ||l2(o,t* -v) 

< L's .11911*2(0,^. ; [1.2 (fi)]3, > || J _ g||L2(0,T- ;[L2(n)]3) • (27) 


4.2 Proof of Theorems Theorems [T}1^ 

Existence. Consider any g > 0. Then, Lemmas |4.If |4.4| imply that operator SET, 

SFT;e, ^ L 2 ( 0 ,r*;E) (28) 


is well defined, continuous and compact for any choice of T* satisfying (20). 


We now claim that, by suitable choices of q and further refinement (i.e., possible reducing) of 
the value of T*, we can prove that SFT(Sq) C Bq. 

Let us choose q > 2 Ls||uo|1l2 , where Ls is from (25). Then, taking into account ([20|, choose 


0 < T* < min 1 , T; , T) 


/ meas 


K.,V 7 


{sm \ 


\ cn c'^ ICsq J 


where > 0 is small enough to ensure that 


LsLFrK,„(r*_„) < 


q 

2 ’ 


(29) 


see Lemmas 


4.3||4.4 for parameters Lf, ^k,v{T(^v)- Then, making use of (25) and (23), we obtain: 


-r .t^F 


Applying Schauder’s fixed point theorem completes the proof of existence in Theorems and 

H o 

Proof of uniqueness. We argue by contradiction. Let (u^^/zi^i) and be two 

distinct solutions to a© for the same initial condition an d for cing term, both defined on the 

this means that and 


4.1 


same time-interval [0,r]. Due to linearity of (|4|) and Lemma 

must be distinct on [0,T]. Hence, since solutions to (|^ are continuous in time with values in 
H, there exists some interval, say, {T,T + t),T > 0,r > 0, such that these two solutions are 
distance at every point of it (in the above sense of continuity) but •) = ■) = u{T). 

Setting (similar to the proof of existence): 


q := max < 


,(i)| 


L2(T,T(i) -V) 1 


,( 2 ) 


Il2(T,T(i) -V) ;2Ls||u^^^(r)||L2| , 
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we choose r* G — T] so that 


0 < T* < min 


/ meas(S'(0)) \ 
\2cqc'^ICs q) J 


with T*^ as in (29) with q/2 in place of q/2. In particular, by defining the operators S,F,T 
as in ([l7|)-([I^ with the time interval [0, T*] replaced by [T, T + t*] and with Bq replaced by 
Bq C L^T, T + T* ]V), we can repeat the previous calculations yielding SFTyB, C Bq and that 
Lemmas |4.2tj4.4| apply. 


Since both and are fixed points of SFT on Bq, we have from (26) in the case when 
d = 2 (and we use (27) for d = 3): 


- ■u(^)||L2(T,T+r 


. ;V) = ||SFTu(l) - SFTu(2)||L2(T,T+r* ;V) 




(T.T+.. ||FTu( 1) - FTu(2)|lL2(T,T+r* ;[L2(n)]^) 


In turn, from (24) it follows that 

||FTu(i)-FTu(2)|lL2(^,r+,.,[L2].) 

Finally, (22) yields: 




< i'x 


cn Ks q v7* 


I - 


- U*^^^||l2(t,T+t* -V)- 


meas{S'(0)) 


Combining the above chain of estimates allows us to conclude that 

— M^^^||l2(T,T+t* ■,V) < ~ W^^^||L2(T,T+r* ;V) ) 

where T(t*) —>■ 0 as t* tends to zero. Thus, and must coincide on [T,T + r»] for some 
r* G (0,T*). Contradiction, o 


5 Proofs of Lemmas 14.11 - 14.4 

Below we use the symbol |a|, a G to denote ||a||Rd. 


5.1 Proof of Lemma 14.11 

Fix T > 0 and u G L'^{0,T; [iJQ(n)]'^) and assume ([^. 

Step 1. Let ho and ICs as in (H2), and choose Tq such that 


0 < To < min < T, 


meas(5'(0)) 


meas(S'(0)) ( 

4 II'^IIl2(o,t ;[L2(n)]2) y 2 Cn Cq ICs II'U'||l2(o,t ;[i/i(a)]2) ^ 


where cn is the constant in (51) and c(^ is the constant in (52). Denoting 

^ ho 


Bho/2 ■— 1 C £ C([0,To] ;M ) ; ||CIIc([o,To]; 


( 30 ) 
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we define, for each index * = 1,..., a function Di: Zi^o + ^?io /2 —>■ C([0, Tq] ; R‘^) by setting 
for all w G + Bh ^/2 


Di{w){t) := z^^o 


meas(S'(0)) Jq Js{w(t)) 


i{T,x)dxdT, tG[0,To]. 


Step 2. We claim that each Di takes values in Zi^ + 5/io/2 and is a contraction mapping. 
Proof of claim. For all t € [0, Tq], we have, : 


|AM(t) - z^,o| < —JtvTva I / / u{T,x)dxdT\ 

measf*(,Ujj Jq Js(w(t)) 


< 


meas(S'(0)) Jq Js(wir)) 

V% I, II ^ ho 

/ li^llL^(0.T;[L2(n)]<i) < ^ , 

A/meas(o(0)) ^ 


where we have used (HI), Schwartz inequality and ( [M| . This proves that Di'. Zi^ + i?/io /2 
Zifi + Bha/2- 

Moreover, if we take wi,W 2 € Zi^ + then for all t G [0, Tq]: 


\Di{wi){t) - Di{w2)it)\ < 


1 


meas(S'(0)) 

1 

meas(S'(0)) 


n u{T, x)dx — / u(r, a;) dx dr 

’{wi{t)) J S(w9.(t)) 


S(«I2(t)) 


n u{T, x) (^i.T(a;) - 6.r (a;)) dx dr 
! 


where we have denoted with the characteristic functions of S(wi{T)). 

For all T € [0,t] such that wi{t) ^ W2{t), by setting /i(r) := wi{t) — W2{t) and denoting with 
i/(r) := r]{h{T)) the unit vector given by (H2) for /i(r) G Bh.g{0) \ {0}, we can consider the 
sections of u{T,-){fi^r — ■Ca.r) corresponding to y G ^u{t) is the projection of fl on the 

hyperplane orthogonal to rj). For these sections, there holds, see notations 


[ u{t,x) {fi^r{x) - f 2 ,r{x)) dx = f ( f a) (6.t “ («) da ) ^{y), 

JQ. J 

where T-D~^ stands for the (n — l)-dimensional Hausdorff measure over Borel sets of K". Then, 
due to (511, for each component Uj of u = (ati,..., Ud): 

/ Uj{t,x) - £,2,t{^)) hx 

an 

ll(uj)^(^)(r,-)l!L“(n“(^,) (6,r - 6,r)5;(^) (a) da^ dn‘^-\y) 

(Ci.r - (a) dc^dn'^-\y). 


< 


< CQ 


Now notice that, (fi^r — £.2,t) coincides with the characteristic function of the symmetric differ¬ 
ence S{wi{t)) a S{w 2 {t)), and that we have 

S{wi{t)) AS{w2{t)) = {wi{t) + S{0)) A{w2iT) + S{0)) = {h{T) + S{0)) AS{0), 
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thanks to (HI). Thus, the above integral over is actually the measure of the section along 
r]{h{T)) of this symmetric difference and, owing to (H2), we obtain: 



(a) 


da < JCs ki(r) - W 2 {t)\ < JCs ||wi 


^2|!c([O.To] ;R-^) ■ 


(31) 


Observing that, for times r G [0,t] such that wi(t) = W 2 {t), the symmetric difference is empty 
and — ^ 2 ,T = 0 and, making use of (521, we can deduce next that 


\D^{wi){t) - A(w 2 )(i)P < 


1 


(meas(S'(0)))2 

d „t 


n u{T, x) {^i,t{x) - ^ 2 ,t{x)) dx dr 
! 


(meas(S'(0)))2 ^ 


X! / / {^^,r{x) - ^ 2 ,t{x)) dx dr 

Jo Jn 


< 


(meas(S'(0)))^ 


E{/7 

j — 1 ^ ^ 

^ 7 “ ^2,r)E) (“) dn‘^~^{y) dx dry 


^ r Cg l^s ll^^l ~ '»^2|ic([0,To] ;R-^) j 2 


meas(S'(0)) 


i=i 


0 Jo, 


(t) 


dn {y) dr 


^Cn/Csllwi -W2 ||c([0,To];ROt2 r° / II ! Ml M 

^ i_ _ -1 E / cn||uj(T,-)||m(n)C^T 

,=1 Ko 


< { 


meas(S'(0)) 

Cn Cj.; II Wl - W 2 | 1 c([ 0 .To] ;Ii 


meas(S'(0)) 


1=1 


■) 11 ^ 1 ( 0 ) 


dr 


< { 


Cn Cq Ailg II Wi — W 2 ||c([ 0 ,To] ;1 
meas(S'(0)) 


-Y Tq, ||u||l2(o^To 


(32) 


Therefore, taking the supremum over t G [OjTq] and recalling (30), we conclude that 


11 ^( 1111 ) - £>j(w2)||c([0,To] IR'*) < 2 11^1 “'*^2||c([0,To] ;R-i) J 
and the proof of the above claim is completed. 

Step 3. By applying the contraction mapping theorem to each Di, we obtain uniquely deter¬ 
mined fixed points z^,... ,z’^ in Zi^ + S/to /2 C C([0,To] ;K‘^). Since the choice of Tq in (30) 
is time invariant, we can repeat the same argument for intervals [nTg, (n -j- l)To], for n > 1, 
until we reach T. Thus, we obtain curves z^,... ^z’^ in C([0, T] ;K'^) which satisfy the integral 
representations 


A A) — ^i,0 + 


meas(S'(0)) Jq Jsiz-ir)) 


i{t, x) dx dr , 


t e [0,T], i = , 
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and, hence, they are solutions to Q. By continuity in time of z*(t)’s, it follows from Q that 
we can choose T* G (0,T] small enough to have 


\z:it)-z*it)\>2r 

for all required indeces and for alH € [0,T*], so that ([^ is satisfied. 

Step 4. Recalling that ||u||l 2 (o,To ;[L 2 (n)]rf) < C'n||'a||L2 (o.Tn due to Poincare inequality, it 

is immediate to see that the last statement in Lemma |4.1| 

meas(5'(0)) hg 


0 < Tfc < min < T, 




follows by replacing ( |30[ ) with 
2 ' 


meas 


{sm \ 


2cnc'^ICskJ 


(33) 


This concludes the proof, o 

5.2 Proof of Lemma 13.11 

Let us start from the term fei- It is immediate to see that, for any fixed {t,x) G (0,T) x Q, 
there holds 


N 


\Zi{t) - Zi-i{t)\ - ii-i 

— 2 max {||zi — Zi_i|lc([ 0 ,T] ;K‘i) + ^i-l} I ^ l'^j-l(OI ) • 


Thus, the first inequality in Lemma |3.1| follows. 
Similarly, for frot, 2 d we have 


N-l 




i=2 


N-l 




2 = 2 


< 2 _^^max^^ < \\zi — 2i-i||c([o,T] + 


\Zi - Z^-l\ 


C([0.T] ;B‘i) 


^Af-1 


2r 




because x Ax is an isometry. Hence, for every t G (0, T) we deduce 
y I/rot, 2d(i,a:)pda:^ < 2 v^meas(H) 


X ^ max < \\Zi — ^i-i||c([o,T] iB-i) + 

2 — . .,i\ 1 I 


\Zi - Zi-i\ 


C([0,r] ;I 


2r 


so that the second inequality in Lemma [3. 1 1 follows as well. 

It remains to prove the third inequality about frot, 3 d- H is immediate to observe that for 
i = 2,..., TV — 1 there hold: 


Pi[t]{zi-i{t)-Zi{t)) = ^i{t)-\zi-i{t)-Zi{t)\, Qi[t]{zi+i{t)-Zi{t)) = ^^{t)■\z^+l{t)-Zi{t)\, 
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where ^i{t) := — Zi{t)) x {zi+i{t) — Zi{t))\. Thus, we can proceed as in the 2-D case to 

obtain 


Af-l 




+ +Ci+i(i,a;)) 


^AT-l 


< 4 [ ^ ki-i(i)l ) ~ '^*-1 


1C([0.T] ;R‘i) > 


where we have used the fact that ^i{t) < \zi{t) — Zi-i{t) \ \zi{t) — Zi+i{t)\. Therefore, the required 
inequality follows as in the case of frot, 2 d- This completes the proof, o 


5.3 Proof of Lemma 14.21 


The operator T is well defined for all times T* > 0 satisfying the assumptions of this lemma, 
thanks to Lemma 4.1 and its images satisfy the properties ([^. Therefore, it remains to prove 
that it is bounded, continuous and compact. 

Step 1: Proof of the continuity of T. 


The proof of this part is similar to that of Lemma 4.1 Let G Bq be fixed. Denote 

= z^^'^ = {zi \ ..., z^jp), j = 1, 2. Then we have for all t G [0, T*]: 


zi^\t)-zr{t)\< 


( 2 )/ 


1 


< 


meas(S'(0)) 

1 

meas(S'(0)) 


n u^^\T,x)dx— [ x) dx dr 


1 


/O 
ct 

■ 

/o JQ. 

rt 


x) — x) ^i^^{x)dx dr 


meas(S'(0)) 


[ f u‘'^\t,x) {^i^r{x) - ^2 ,t{x)) dx dr 

Jo Jq 


where we have denoted with ■Cj,r(a;) the characteristic functions of S{zI‘^\t)). By applying 
Holder inequality to the first term and by repeating for the second term the argument used to 
derive (32)-(331 in Lemma 4.1 we obtain: 




( 2 )/ 


meas(S'(0)) 




-b 


Cfi Cq K-s 'Jt* 
meas(5'(0)) 


,( 2 ) 


(1) J2)| 


M l|L 2 (o.T*;y)lDi - 2 :^ ||c([o,r*] 


< Cn 




+ 


meas(S'(0)) 
cncJQlCsqVT* (i) 


- u(^^|lL2(o,r*;y) 


z' — z' 


(2) I 


|C([0.T*] ;R‘i)- 


meas(S'(0)) 

After taking the supremum over [0, T*] and rearranging the terms, we arrive next at the estimate 
cqc^qICs qVJ^\ II (1) 


1 - 


~ 4^^IIc([ 0,T*] ;R^) < C'n ' 


meas(S'(0)) 


I T* 
meas(«S'(0)) 
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implying (221 and the continuity of T, due to (20). 


Step 2: Proof of the boundedness of T. 

Let u be fixed in Bq and consider any i = 1,..., N. If we denote with Zi the i-th component of 
Tit, we obtain for all t G [0, T*] that 


\Ziit)\ < \Zifi\ 


1 


meas(S'(0)) 


Thus, 


0 JS(Zi{T)) 


u{t, x) dx dr 


Tm||[c([o,t*] iR'*)]" < j l^^ol + C'n 


< ko| + 




meas(5'(0)) 




meas(S'(0)) 


l^l!L2(0.r* -H)- 


T l|w|!L2(0.r* ;V) 


Step 3: Proof of the compactness of T. 

We will show that any bounded sequence from Bq is mapped by T to a sequence 

{zU) which has a converging subsequence. To this end, we will prove that 

is bounded and equicontinuous. Then Ascoli-Arzela theorem will provide the de¬ 
sirable result. 

The uniform boundedness of follows immediately from Step 2. To prove equiconti- 

nuity, Hx t,t + h G [0, T*] and assume that h > 0 (the case h < 0 is analogous). Then, for every 
index i = 1,..., A^, the i-th component of satisfies: 


z^^\t + h)-zl^>it)\ = 


U)( 


1 

meas(S'(0)) 




u^^\t,x) dxdr 


< 


A/meas(S'(0)) A/meas(S'(0)) 


This completes the proof of Lemma 4.2 o 


Remark 5.1 The results similar to Lemmas 4-1 4-^ were proven in m-ina and m 

for swimming models similar to but for the case of nonstationary Stokes fluid equations. 

However, m-m and m deal with the case when u is more regular, namely, when u G 
L2(0,T;[L°°(Oy) C L2(o,T;[ij2(II)]d). 


5.4 Proof of Lemma 14.31 


The operator F is well defined under the assumptions of Lemma [4~3l thanks to Lemma 4.2 


In turn, the estimates (23)-(24| can be derived by straightforward transformations of explicit 
elementary expressions in the forcing term (10) (see also the proof of Lemma |3.1[ ). Moreover, 
the results very similar to Lemma 4.3 were derived in all technical detail in |16|. [T7] (Section 
12.4) for the 2-D case and in (TH] for the 2>-D case for swimming models similar to (|^-([^ but for 
the case of nonstationary linear Stokes fluid equations. The difference between the forcing terms 
in |16j-|17j and |19j and in this paper is that in the former nfs were positive numbers and Vi’s 
were elements of L°°{0,T), while in this paper Kfs and vfs lie in Lf{{),T). Furthermore, in the 
case of Z-D swimming models the formula for rotational forces in m has a more general shape 
but designed to work only when the respective triplets {zi-i, Zi, Zi+i} are not in the aligned 


position (see subsection 3.2 in the above). 
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5.5 Proof of Lemma 14.41 

The classical results in Theorems and about solutions to incompressible Navier-Stokes 
equations cited in Section 6 below, ensure that operator S is well defined. Let us also 
remind the reader that the initial datum uq is fixed in our proofs. 


Proof of the boundedness of S. 

For any / € L^(0,T* ; [L^(r2)]'^), estimate (25) follows immediately from (47) in the 2-D case 
and from (|4^ in the 3-D case. 


Proof of the continuity of S. 

Step 1. Consider first the case when d = 2. Let (ui = Sfi,U 2 = S/ 2 ) be a pair of solutions 
to uncoupled § generated, given uq, by some /i ,/2 G L^(0,T* ; [L^(f2)]^). To derive ( [^ , it is 
sufficient to establish the following estimate: 


||U1 - U2 


2 

L2(0.T* -V) 


< 




!C{||U1| 


L^(0,T* ;V)’" 


"L^(0,T* ;V) 


Wfi 


/2II 


2 

L2(0,T* ;[L2(a)]2) • 


(34) 


where Cq is the constant from Poincare inequality. 

Step 2. From the proof of Theorem III.3.2 in [30], we obtain that for any t € (0,T*) the 
difference Ui{t) — U 2 {t) satisfies the following equality (see in particular (3.65) in (SO] Chapter III] 
or (5), page 145 in g^): 


where (•, •) 


^ Q WMt) - Mt)\\h'^ + V ||ui(t) - M2(t)|i^i 

= (/l(i) - /2(t),Ml(t) - W2(t))L2 

- b{ui{t) - U 2 {t),U 2 {t),Ui{t) 

l 2 stands for the scalar product in L^(r2). Hence, we have 
^ (||ui(t) - U 2 {t)\\l 2 ) -k 2v ||ui(t) - U 2 (t)||^l 


U2it)) , 


(35) 


< 2||/i(t) - /2(t)|lL2 \\uiit) - U2{t)\\ 


L2 


b{ui{t) - U 2 {t),U 2 {t),Ui[t) - U 2 {t)) 


< ‘2.Cn\\fl{t) - f2{t)\\-L2 \\ui{t) -U2{t)\\Hi 
+ 2V2 ||u 2 (f)||m \\ui{t) - U 2 {t)\\i ^2 ||ui(t) -U 2 (t)|| 




< 


\\fl{t) - f2it)\\l2 + ^ ||ui(t) - U2(t)f 




-.\\Mt)\\m WMt) - Mt)\\h + WMt) - Mt)\\m 


m ’ 


where Cq is the constant in Poncare inequality and, to evaluate |5(-,-,-)|, we used (46) from 
Section 6. Therefore, there holds: 


— (||Ml(t) - U2(t)|!L2) + y ||ui(t) - U2{t)\\ 


m 


< 


2C^ 


l/lW-/2(t)l 


2 

L2 


\u 2 ml. \\ui{t)-u2it)\\l2. (36) 
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Owing to (36), we immediately see that 


Hence, an application of Gronwall’s inequality yields that 

\\ui{t) - U2{t)\\l2 < ^ Wflis)-f2{s)\\l2ds-exp(^^\\u2\\l2(^o^t.y^^ , tG[ 0 ,T*]. ( 37 ) 

In turn, integrating (36), we deduce from (37) that 

v\\ui — U2\\i2^Q rp, .y^ < ||ui(T ) — U2{T )||l2 + v\\ui — M2|Il2(o,t* -v) 


< 


2C? 


< 


f 4 

/ ll/i(s)-/2(s)||Lds+ max ||mi(s)-M2(s)| 1L • - ||M2|lL(0,T*;y) 

Jo s6[0.T*] ly 

^ J 11/1(5) - / 2 (s)| 1 l 2 ^ lk 2 |lL 2 ( 0 ,T* ;V) ' 


From here, making use of estimate (1 + xe“) < for a: > 0, we arrive at 


l“i - W2 ||l2(o,t* -y) - — 2 ^ •!'> Il/i - /2|lL2(o,r* 


[L2(n)]0) 


(38) 


implying (26). This completes the proof in the 2d case. 


Remark 5.2 An estimate similar to (34) can also be found in 121}/ . see Theorem 11 and esti¬ 
mates ( 45 ) and ( 48 ) on pp. 170-171. 

Step 3: Case d = 3. Similar to Step 2, consider any two solutions Ui = Sfi,i = 1,2 to ([^, 
generated by some /i ,/2 in L^(0,T* ]H). Combining estimate (11) from [^, p. 147, namely: 

llV’IU^'Cn) < '^'0 G ^o(^): where D* is some positive constant, 

with the fact that U 2 G C'([0,T*]; V) (see Theorembelow) allows us to define 


C2= max |M2(t)|L4 < 3^/'‘£)*||M2||c([o.T*];y), 

tG[0,T*] 'in VI . i. ; 


(39) 


where 


U2(t)|L4 



U2 = (u21,M22,M23)- 


To evaluate b(^ui{t) — U 2 {t),U 2 {t),ui{t) — U 2 {t)) in the 3-D-case, we will use the chain of estimates 
from HD, p. 145 as follows. 


b{ui{t) - U2{t),U2{t),Ui{t) - U2it)) 


< V3|m2(0|l4 hl(t) -U2{t)\\H/ - U2it)\-Li 

< \f 3 C 2 \\ui{t) - U 2 {t)\\H^ I Ui{t) - M2(i)|L4- 
Then, making use of estimate (12) from [^, p. 147 and (5) on p. 10, namely: 

o3/4 

ll/’llL'‘(n) < ^llV'llffi(n) + -^ll/’llL2(n): 4’ G <5 > 0, 
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we obtain as in [^, p. 145: 


b{ui{t) - U 2 {t),U 2 {t),Ui{t) - U 2 {t)) 


< V3C2\\ui{t) - U2{t)\\ 


m 


e\\ui{t) - U 2 {t)\\Hi + Ce\\ui{t) - U2(t)||[L2(Q)]3 




< 2V3C2e\\ui{t)-U2it)\\Hi+ViC2^\\ui{t)-u2{t)\\fL2^Q)]3, e>0, 

where = 3(3)^/"‘/e^. 

Select now e = (2v^C2})“^(j 2/4), then: 


, 


b{ui{t) - U 2 {t),U 2 {t),Ui{t) - U 2 {t)) 


< -limit) - U2(t)||^. + -^llmit) - r.2(t)||f,.2(0)]3, 


4 “0 jy 

|2 


< j\\ui(t) - Mt)l\Hi + Delimit) - U2(t)||[L2(o)]3, 


where we can set, due to (50) and (39) 


(40) 


£>* = CI/I (||uo||y + ||/2 ||l2(0.Ti ;[L2(n)]3))® , 


(41) 


and C depends on v. Similar to (35), we obtain: 


I (||ui(t)-U2(t)||^2)+2H|ui(t)-U2(t)||^i 


(42) 


< 2||/i(t) - /2(t)|lL3 ||ui(t) - M2(t)|lL2 +2 b{uiit) - U2it),U2it),Ui{t) - U2it)) 


< 


2C, 


Whit)-f 2 mi.+-limit)-U 2 m 


-||(ui(t) - U2(t))||il + 2D*||(Mi(t) - U2(t))||[i2(0)]3 . 


Respectively, instead of (36), we have: 


— (||Ml(t) - U2(t)|lL2) + ||ui(t) - U2(t)||ffi 


< 




and, furthermore, in place of (38): 



u 

38 

I: 


Whit) — /2(0liL2 + 2Il*|j(Mi(<) — U2(t))||p^2(Q)]3 , (43) 


4 (^2 

|mi — U2||L2(o_r. ;y) < 2 ® * ll/l ~ /2 ||l 2(0.T* ;[L2(n)]2)) : 


(44) 


implying (27) in view of (411. This completes the proof of Lemma 4.4 o 
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6 Appendix A: Solutions to uncoupled Navier-Stokes equa¬ 
tions 

In this section we will remind the reader several classical results relevant to the well-posedness 
of Navier-Stokes equations which are used in the proofs of the main results of this paper. 
Assume that in ([^: 

f€L^{0,T;[L^{D)f), uq G [L2(U)]". 
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Following [301 Chapter III], we introduce the following trilinear continuous map on [iJQ(n)]‘^ x 
[Hq^H)]^ X [iJQ(n) (see Lemma II.l.l in [3D], page 161-162 for d = 2,3, page 159): 


b{u, V, w 



Ui (DiUj) Wj dx. 


(45) 


Furthermore, due to Lemma 3.4 (page 292) in (SD] Chapter III], for d = 2, for all u,v,w G 
[ido(n)]^ there holds 

|6(w,z;,'u;)| < ^2 (UmUj^i 1 |w 1 |l 0 " (Iklli/i IklkO' • (46) 


Our results concerning the well-posedness of the swimming model will heavily rely on the 
following classical results about existence and uniqueness of weak solutions for ([^. 


Theorem 4 (Theorems III.3.1 III.3.2 in |30) . pages 282, 294) Let T be any positive 
number, d = 2, / G L^(0, T ; [L^(f2)]^) and uq G H. Then, there exists a unique solution 
wG C([0,r];7d)nL2(0,T;l/) of (§ on [0,r] and the following estimate holds (see estimates 
( 45 ) and ( 48 ) in Theorem 11 \21^ . pp. 170-171 as well as our estimates (36)-(38) when ui = u 
and U 2 , /2 = 0): 


l|w||c(0.T;ff) + I1 w|1l2(0,T;V) < Ls { ||wo| 1 l 2 + ||/||L2(0,T;[L2(n)]2)) 


(47) 


for some positive constant Ls ■ 


Once solution u is found, the following result allows us to retrieve p. 

Proposition 6.1 (Theorem V.1.7.1 in [27| . page 295) Under the assumptions of Theorem 
there exists a function V G L^(0, T ; L^(0)) such that p = Vt is an associated pressure of u 
and 

Vp = f — Ut vAu — (u ■ V)u 
in the sense of distributions in (0,T) x LI. 


Theorem 5 (Theorem 9, [20] . p. 203 and Lemma 9, p. 194) Let T > 0 be any given 
positive number, d = 3, f G L^(0,r; [L^(0)]^), uq G V and dLl be of class C^. Then for any 
Ti satisfying 


TiG( 0, iG), K —mini T,C° ||uo||^i-I-||/||L2(o_r.[L2(f2)]3) 


(48) 


(C° depends on LI and v only), ([^ admits a unique solution {u,p) on (0, Ti) such that ut, Am, Vp G 
[L^((5ti)]^ and u G C([0,Ti];y). Furthermore, the following estimates holds (see (42)-(44) 
when Ml = M and U 2 , f 2 = Oj; 


l|w||c([0,Ti];ff) + ll'u||L2(0,Ti;y) < (I|uo||l2 + II/IIl2(0,Ti ;[L2(n)]3)) 

for some positive constant Lg > 0 and (see Lemma 9, p. 194, (55)): 

l|Ml|c([0.Ti];y) < is (llMollffi + ||/||l2(0,Ti ;[L2(n)]3)) 


(49) 

(50) 


for some positive constant Ls > 0. 

The cited Theorem 9 in |20|, p. 203 admits the following modification. 

Theorem 6 (Additonal regularity, |21| . Theorem 17, pp. 183-184) //mq G [7d^(f2)]^ f) T, 
then in Theorem^solution u lies in [id^d(Qj,^)]3 Q(27([0, Ti]; C). In turn, ifuo G [H^{Ll)Y 0V, 
then u G C\C{[Q,Ti\-,V) Theorem^ 
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7 Appendix B: Sections 


In the proof of Lemma [O] in the above we have used the following results (see assumption (H2) 
for notations). 


Proposition 7.1 (Proposition 3.105 in [2]) Let p G [l,cx)), w G L^(fl) and v be a unit 
vector in K". If we define the generalized directional derivative of w in the (unit) direction v 
as a distribution G ll'(fl) such that 


4>T)^wdx = — [ w ^ dx 

! Jn diz 


yf>GCT{n), 


and assume that T)i,w G then for every y G we have G and Dw^ = 

(D^w)l^, i.e., the differential of the section w^, is the section of the directional derivative 
In particular, if w G then wj) G for all vectors v, \v\ = 1 and all y G Lly. 


In the proof of Lemma 4.1 we use the following corollary, where we denote by "H" ^ the (n — I)- 


dimensional Hausdorff measure over Borel sets of . 


Corollary 7.1 Let w G IlQ{n). Then, the following properties hold. 

(i) There exists a constant cq > 0 such that, for all vectors v, \v\ = 1 and all y G 

lk^llL“(n«) < cn \\wI\\h^(qi) ■ 

(ii) There exists a constant Cq > 0 such that, for all vectors v, \v\ = 1, we have: 

\\wl\\H^{Ql)dn'^-\y) < c'n ||w||ffi(n). 


we have: 

(51) 

(52) 


Proof: Part (i) is an immediate consequence of Sobolev embedding theorems for bounded sets 
C K, see e.g. [1], and of the equivalence between the norm we defined on and the 

usual one on Indeed, the Sobolev constants depends only on the measure of the sections 

and the latter is bounded by the diameter of Q, independently on the direction v and on 
the point y G it l,. 

Concerning (ii), it is enough to apply Holder’s inequality and recall Proposition 7.1 to obtain: 

1/2 


= v/H"-i(^^.) 


'a/ 


\T)wl{a)G da ] dTff {y) 


[ ([ da) dn-\y) 

Jq.„ \Jni J 


1/2 


-I 1/2 


\Dvw{x)\^ dx 


Now, notice that the 1st term in the resulting expression, can be estimated, 

independently on the chosen direction v, by the (n — l)-dimensional volume of a ball of radius 
i? > 0 with R containing Cl. For the second term we have, in turn: 


ID, 


n p 

;(x)p dx < ^ / \Diw(x)\^ dx = 

dn 




whence the desirable conclusion follows, o 
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